A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar, and an active nucleus. We present a detailed study of the spatial variation of the far-infrared (FIR) [C II]158 μm and [O I]63 μm lines and mid-infrared H_2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star-forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([C II]158 μm+[O I]63 μm)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7 μm PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [C II]158 μm/PAH(5.5-14 μm) is found. PAHs in the ring are responsible for a factor of two more [C II]158 μm and [O I]63 μm emission per unit mass than PAHs in the Enuc S. spectral energy distribution (SED) modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high-intensity photodissociation regions (PDRs), in which case G_0 ~ 10^(2.3) and n_H ~ 10^(3.5) cm^(–3) in the ring. For these values of G_0 and n_H, PDR models cannot reproduce the observed H2 emission. Much of the H2 emission in the starburst ring could come from warm regions in the diffuse interstellar medium that are heated by turbulent dissipation or shocks

H2 formation and excitation in the Stephan's Quintet galaxy-wide collision

Context. The Spitzer Space Telescope has detected a powerful (L(H2)~10^41 erg s-1) mid-infrared H2 emission towards the galaxy-wide collision in the Stephan's Quintet (SQ) galaxy group. This discovery was followed by the detection of more distant H2-luminous extragalactic sources, with almost no spectroscopic signatures of star formation. These observations set molecular gas in a new context where one has to describe its role as a cooling agent of energetic phases of galaxy evolution. Aims. The SQ postshock medium is observed to be multiphase, with H2 gas coexisting with a hot (~ 5 10^6 K), X-ray emitting plasma. The surface brightness of H2 lines exceeds that of the X-rays and the 0-0 S(1) H2 linewidth is ~ 900 km s-1, of the same order of the collision velocity. These observations raise three questions we propose to answer: (i) Why H2 is present in the postshock gas ? (ii) How can we account for the H2 excitation ? (iii) Why H2 is a dominant coolant ? Methods. We consider the collision of two flows of multiphase dusty gas. Our model quantifies the gas cooling, dust destruction, H2 formation and excitation in the postshock medium. Results. (i) The shock velocity, the post-shock temperature and the gas cooling timescale depend on the preshock gas density. The collision velocity is the shock velocity in the low density volume filling intercloud gas. This produces a ~ 5 10^6 K, dust-free, X-ray emitting plasma. The shock velocity is smaller in clouds. We show that gas heated to temperatures less than 10^6 K cools, keeps its dust content and becomes H2 within the SQ collision age (~ 5 10^6 years). (ii) Since the bulk kinetic energy of the H2 gas is the dominant energy reservoir, we consider that the H2 emission is powered by the dissipation of kinetic turbulent energy. (Abridged)Comment: 19 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics Minor editing and typo

H_2 formation and excitation in the Stephan's Quintet galaxy-wide collision

Context. The Spitzer Space Telescope has detected a powerful (L_(H_2) ~ 10^(41) erg s^(-1)) mid-infrared H_2 emission towards the galaxy-wide collision in the Stephan's Quintet (henceforth SQ) galaxy group. This discovery was followed by the detection of more distant H_2-luminous extragalactic sources, with almost no spectroscopic signatures of star formation. These observations place molecular gas in a new context where one has to describe its role as a cooling agent of energetic phases of galaxy evolution. Aims. The SQ postshock medium is observed to be multiphase, with H_2 gas coexisting with a hot (~5 × 10^6 K), X-ray emitting plasma. The surface brightness of H_2 lines exceeds that of the X-rays and the 0-0 S(1)H_2 linewidth is ~900 km  s^(-1), of the order of the collision velocity. These observations raise three questions we propose to answer: (i) why is H_2 present in the postshock gas? (ii) How can we account for the H_2 excitation? (iii) Why is H_2 a dominant coolant? Methods. We consider the collision of two flows of multiphase dusty gas. Our model quantifies the gas cooling, dust destruction, H_2 formation and excitation in the postshock medium. Results. (i) The shock velocity, the post-shock temperature and the gas cooling timescale depend on the preshock gas density. The collision velocity is the shock velocity in the low density volume-filling intercloud gas. This produces a ~5 × 10^6 K, dust-free, X-ray emitting plasma. The shock velocity is lower in clouds. We show that gas heated to temperatures of less than 10^6 K cools, keeps its dust content and becomes H_2 within the SQ collision age (~5 × 10^6 years). (ii) Since the bulk kinetic energy of the H_2 gas is the dominant energy reservoir, we consider that the H_2 emission is powered by the dissipation of kinetic turbulent energy. We model this dissipation with non-dissociative MHD shocks and show that the H_2 excitation can be reproduced by a combination of low velocities shocks (5-20 km s^(-1)) within dense (n_H > 10^3 cm^(-3)) H_2 gas. (iii) An efficient transfer of the bulk kinetic energy to turbulent motion of much lower velocities within molecular gas is required to make H_2 a dominant coolant of the postshock gas. We argue that this transfer is mediated by the dynamic interaction between gas phases and the thermal instability of the cooling gas. We quantify the mass and energy cycling between gas phases required to balance the dissipation of energy through the H_2 emission lines. Conclusions. This study provides a physical framework to interpret H_2 emission from H_2-luminous galaxies. It highlights the role that H_2 formation and cooling play in dissipating mechanical energy released in galaxy collisions. This physical framework is of general relevance for the interpretation of observational signatures, in particular H_2 emission, of mechanical energy dissipation in multiphase gas

High Performance Morphological Filtering of Cirrus Emission from Infrared Images

The project was designed to explore new morphological filtering, techniques for the removal of foreground Galactic 'Cirrus' emission from NASA Infrared Astronomical Satellite (IRAS) data, especially at 100 microns, using parallel processors as the main engine for achieving this result. The ultimate aim was to provide NASA with completely filtered data by the end of the grant period for the entire IRAS 100 and possibly 60 micron database. If successful, the filtered data would reveal many new sources of IR emission, especially at low galactic latitudes, which had previously been heavily confused with diffuse Galactic 'cirrus'

Cloud fluid models of gas dynamics and star formation in galaxies

The large dynamic range of star formation in galaxies, and the apparently complex environmental influences involved in triggering or suppressing star formation, challenges the understanding. The key to this understanding may be the detailed study of simple physical models for the dominant nonlinear interactions in interstellar cloud systems. One such model is described, a generalized Oort model cloud fluid, and two simple applications of it are explored. The first of these is the relaxation of an isolated volume of cloud fluid following a disturbance. Though very idealized, this closed box study suggests a physical mechanism for starbursts, which is based on the approximate commensurability of massive cloud lifetimes and cloud collisional growth times. The second application is to the modeling of colliding ring galaxies. In this case, the driving processes operating on a dynamical timescale interact with the local cloud processes operating on the above timescale. The results is a variety of interesting nonequilibrium behaviors, including spatial variations of star formation that do not depend monotonically on gas density

Multi-waveband observations of colliding galaxies

Colliding galaxies represent a major challenge to both theorists and observers because of the large variety of phenomena which are expected to come into play during the interaction. Strong gravitational fluctuations may drive non-linear waves and instabilities throughout the stars and gas leading to enhanced star formation, nuclear activity and ultimately a mixing of the morphological components of the original galaxies. One relatively uncomplicated class of colliding galaxy where stellar waves play an important role in star formation are ring galaxies. Ring galaxies are probably formed when a companion galaxy passes through the center of a disk system driving circular waves through the disk (Lynds and Toomre 1976, Toomre 1978, Struck-Marcell 1990). Off-center collisions can generate non-circular waves and can be loosely described as banana-shaped although they may exhibit more complex forms as the waves expand into the disk. The propagation of such stellar and gaseous waves through the disk leads to enhanced star formation (e.g., Appleton and Struck-Marcell 1987a; Jeske 1986) and provides a unique probe of the response of the interstellar medium (ISM) to a propagating wave (see Appleton and Struck-Marcell 1987b). Here, the authors report results for 3 systems; the irregular ring Arp 143 (=VV 117); Wakamatsu's Seyfert ring (A0959-755; see Wakamatsu and Nishida 1987) and the brighter member of the pair of ring galaxies comprising of AM 1358-221. The most complete multi-wavelength data is for Arp 143. Optical charge coupled device (CCD) observations made with the 60 inch Palomar telescope at BV and r band, near-IR images at J (1.25 microns), H (1.65 microns) and k (2.2 microns) bands from the infrared camera (IRCAM) InSb array camera on the 3.8m United Kingdon Infrared Telescope (UKIRT) telescope and very large array (VLA) observations at 20cm in both the neutral hydrogen line and radio continuum are described. The observations of Wakamatsu's ring and AM 1358 were made only in the near-IR, and a comparison is made with available optical plate material

Enhanced dust heating in the bulges of early-type spiral galaxies

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20–30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 μm. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 ± 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23 (±3)% hotter for morphological types earlier than Sc, and only 9 (±3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar

Observations and modeling of the dust emission from the H_2-bright galaxy-wide shock in Stephan's Quintet

Context. Spitzer Space Telescope observations have detected powerful mid-infrared (mid-IR) H_2 rotational line emission from the X-ray emitting large-scale shock (~15 × 35 kpc^2) associated with a galaxy collision in Stephan's Quintet (SQ). Because H_2 forms on dust grains, the presence of H_2 is physically linked to the survival of dust, and we expect some dust emission to originate in the molecular gas. Aims. To test this interpretation, IR observations and dust modeling are used to identify and characterize the thermal dust emission from the shocked molecular gas. Methods. The spatial distribution of the IR emission allows us to isolate the faint PAH and dust continuum emission associated with the molecular gas in the SQ shock. We model the spectral energy distribution (SED) of this emission, and fit it to Spitzer observations. The radiation field is determined with GALEX UV, HST V-band, and ground-based near-IR observations. We consider two limiting cases for the structure of the H_2 gas: it is either diffuse and penetrated by UV radiation, or fragmented into clouds that are optically thick to UV. Results. Faint PAH and dust continuum emission are detected in the SQ shock, outside star-forming regions. The 12/24 μm flux ratio in the shock is remarkably close to that of the diffuse Galactic interstellar medium, leading to a Galactic PAH/VSG abundance ratio. However, the properties of the shock inferred from the PAH emission spectrum differ from those of the Galaxy, which may be indicative of an enhanced fraction of large and neutrals PAHs. In both models (diffuse or clumpy H_2 gas), the IR SED is consistent with the expected emission from dust associated with the warm (> 150 K) H_2 gas, heated by a UV radiation field of intensity comparable to that of the solar neighborhood. This is in agreement with GALEX UV observations that show that the intensity of the radiation field in the shock is GUV = 1.4±0.2 [Habing units]. Conclusions. The presence of PAHs and dust grains in the high-speed (~1000 km s^(-1)) galaxy collision suggests that dust survives. We propose that the dust that survived destruction was in pre-shock gas at densites higher than a few 0.1 cm^(-3), which was not shocked at velocities larger than ~200 km s^(-1). Our model assumes a Galactic dust-to-gas mass ratio and size distribution, and current data do not allow us to identify any significant deviations of the abundances and size distribution of dust grains from those of the Galaxy. Our model calculations show that far-IR Herschel observations will help in constraining the structure of the molecular gas, and the dust size distribution, and thereby to look for signatures of dust processing in the SQ shock